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Abstract. It has been recently shown that the quantum Boltzmann equations may 
be relevant for the leptogenesis scenario. In particular, they lead to a time-dependent 
CP asymmetry which depends upon the previous dynamics of the system. This memory 
effect in the CP asymmetry is particularly important in resonant leptogenesis where the 
asymmetry is generated by the decays of nearly mass-degenerate right-handed neutrinos. 
We study the impact of the nontrivial time evolution of the CP asymmetry in resonant 
leptogenesis, both in the one-flavour case and with flavour effects included. We show 
that significant qualitative and quantitative differences arise with respect to the case in 
which the time dependence of the CP asymmetry is neglected. 



1. Introduction 

Thermal leptogenesis [U, El E] is a simple mechanism to explain the observed baryon 
number asymmetry (per entropy density) of the Universe = (0.87 ± 0.02) x 10^^*^ [4J. A 
lepton asymmetry is dynamically generated and then converted into a baryon asymmetry 
due to {B + L)-violating sphaleron interactions [5l [6] which exist in the Standard Model 
(SM). A simple model in which this mechanism can be implemented consists of the SM 
plus three right-handed (RH) Majorana neutrinos. In thermal leptogenesis [2] the heavy 
RH neutrinos are produced by thermal scatterings after inflation and subsequently decay 
out-of-equilibrium in a lepton number and CP-violating way, thus satisfying Sakharov's 
constraints [71 [6]. RH neutrinos are also a key ingredient in the formulation of the well- 
known "see-saw" (type I) mechanism [8], which explains why neutrinos are massive and 
mix among each other and why they turn out to be much lighter than the other known 
fermions of the SM. For such a reason thermal leptogenesis has been the subject of intense 
research activity in the last few years with the main goal of providing a quantitative 
relation between the light neutrino properties and the final baryon asymmetry. 

Thermal leptogenesis is based on the assumption that RH neutrinos are efficiently 
generated by thermal scatterings during the reheating stage after inflation. In the scenario 
in which the RH neutrinos are hierarchical in mass, successful leptogenesis requires that 
the RH neutrinos are heavier than about 10^ GeV [H]. Hence, the required reheating 
temperature cannot be much lower [21 UHl El [12] ■ In supersymmetric scenarios this may 
be in conflict with the upper bound on the reheating temperature necessary to avoid 
the overproduction of gravitinos during reheating [13|. Indeed, being only gravitationally 
coupled to the SM particles, gravitinos may decay very late jeopardising the successful 
predictions of Big Bang nucleosynthesis. 

In the resonant leptogenesis scenario [TH |T2] this tension may be avoided. If the RH 
neutrinos are nearly degenerate in mass, self-energy contributions to the CP asymmetries 
may be resonantly enhanced, thus making thermal leptogenesis viable at temperatures 
as low as the TeV. Resonant leptogenesis seems to us a natural possibility. In the 
absence of simple predictive models one expects that left-handed and RH neutrinos show 
similar levels of degeneracy. Indeed, within the see-saw mechanism, quasi-degenerate 
light neutrinos are more naturally explained by quasi-degenerate RH neutrinos, rather 
than by an interplay between Yukawa couplings and the masses of the RH neutrinos. The 
quasi- degeneracy in the RH sector can be easily explained by symmetry arguments, e.g. 
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a slightly broken S0{3) symmetry. 

Resonant leptogenesis is the subject of the present paper. More specifically, we will 
investigate the resonant leptogenesis scenario in view of the recent results achieved by 
studying the dynamics of thermal leptogenesis by means of quantum Boltzmann equations 
[16] (for an earlier study, see Ref. [T7]). 

Let us pause here for a moment and summarize why quantum Boltzmann equations 
are relevant in resonant leptogenesis 0- The generation of the baryon asymmetry occurs 
when RH neutrinos are out-of-equilibrium. Therefore, their abundance and the one for the 
lepton asymmetry are determined by Boltzmann equations. In the classical Boltzmann 
equation approach, every scattering in the plasma is independent from the previous one 
and the particle abundances at a given time do not depend upon the previous dynamical 
history of the system. Quantum Boltzmann equations are obtained starting from the non- 
equilibrium quantum field theory based on the Closed Time-Path (CTP) formulation. It 
is a powerful Green's function formulation for describing non-equilibrium phenomena in 
field theory. It allows to obtain a self-consistent set of quantum Boltzmann equations 
for the quantum averages of operators, e.g. the lepton asymmetry operator, evaluated in 
the in-state without specifying the out-state. The quantum Boltzmann equations have an 
obvious interpretation in terms of gain and loss processes. What is unusual, however, is the 
presence of the integral over time in the scattering terms where theta functions ensure that 
the dynamics is causal. The quantum Boltzmann equations are therefore manifestly non- 
Markovian. Only the assumption that the relaxation timescale of the particle asymmetry 
is much longer than the timescale of the non-local kernels leads to a Markovian description. 
A further approximation, i.e. taking the upper limit of the time integral to infinity, 
leads to the familiar classical familiar Boltzmann equation. The physical interpretation 
of the integral over the past history of the system is straightforward: it leads to the 
typical "memory" effects which are observed in quantum transport theory [TS1[T2|. The 
thermalization rate obtained from the quantum transport theory may be substantially 
longer than the one obtained from the classical kinetic theory. 

Furthermore, and more importantly, the CP asymmetry turns out to be a function 
of time, even after taking the Markovian limit. Its value at a given instant depends upon 
the previous history of the system. If the time variation of the CP asymmetry is shorter 
than the relaxation time of the particles abundances, the solutions to the quantum and 

f For more technical details the reader is referred to Ref. [TO] . 
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the classical Boltzmann equations are expected to differ only by terms of the order of the 
ratio of the timescale of the CP asymmetry to the relaxation timescale of the distribution. 
In thermal leptogenesis with hierarchical RH neutrinos this is typically the case. However, 
in the resonant leptogenesis scenario, at least two RH neutrinos A^i and A^2 are almost 
degenerate in mass and the CP asymmetries from the decays of the RH neutrinos are 
resonantly enhanced if the mass difference AM = (M2 ~ Mi) is of the order of the decay 
rates. The typical timescale to build up coherently the CP asymmetry is of the order of 
1 / AM, which can be larger than the timescale for the change of the abundance of the RH 
neutrinos. This tells us that in the case of resonant leptogenesis significant differences are 
expected between the classical and the quantum approach. 

The rest of the paper is organized as follows. In Section [2] we study the impact 
of the time-dependent CP asymmetry on the final lepton asymmetry, in the one-flavour 
approximation. The numerical results are supported by analytical estimates, for both the 
regimes of strong and weak wash-out. The weak wash-out case is analyzed in greater 
detail since this is the regime where the time dependence of the CP asymmetry becomes 
more relevant. The generalization to more than one flavour is given in Section [3l where 
we discuss the possible wash-out regimes. Again, the analytical formulae confirm the 
numerical solutions. Finally, Section H] contains our conclusions. 

2. Resonant leptogenesis revisited: the one-flavour case 

Our starting point is the SM plus three RH neutrinos Na (a = 1,2,3), with Majorana 
masses Mq. The interactions among RH neutrinos, Higgs doublets H, lepton doublets it 
and singlets Cj {i = e, /i, r) are described by the Lagrangian 

^int = Xa^NJiH + Kc^iH' + ]^M^N^N^ + h.c. , (1) 

with summation over repeated indices. The Lagrangian is written in the mass eigenstate 
basis of RH neutrinos and charged leptons. 

For the sake of simplicity, for the rest of the paper we will restrict ourselves to the 
case in which only the two lightest RH neutrinos are quasi-degenerate. Mi ~ M2 -C M3. 
We first address the dynamics of the system in the so-called one-flavour approximation, 
where Boltzmann equations are written for the abundance of the RH neutrinos and for 
the total lepton asymmetry. This approximation is correct only when the interactions 
mediated by charged lepton Yukawa couplings are out of equilibrium. Supposing that 
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leptogenesis takes place at temperatures T ~ Mi ~ M2, the one-flavour approximation 
holds for Ml > 10^^ GeV. We will include flavour effects later on. 

The Boltzmann equations resulting from the OTP formalism, after taking the 
Markovian limit, are given by [H] 



Y'^ = -2e{z)Y^-\Kz^Ki{z)Yc. (2) 

Here Yjq = Y)v^ ~ Yjy^ denotes the number density of the RH neutrinos per entropy 
density which we assume to be roughly equal since we also take the decay rates T^-^ 
and r^Vj roughly equal; primes stand for derivatives with respect to the "time" variable 
z = Mi/T (T being the temperature); the parameter K = Vn^/H^Mi) = T^^/ H{Mi), 
where H{z) is the Hubble rate, controls how much RH neutrinos are out-of-equilibrium; 
Y^ = {l/Ag^)z'^K2{z) is the equilibrium number density of RH neutrinos (being (7* 
the number of relativistic degrees of freedom in the plasma); for consistency, Y^^ has 
been computed using the Boltzmann distribution; Ki and K2 are the modified Bessel 
functions of the first and second kind, respectively. In Eq. ([2]) we have neglected for 
simplicity the contribution of AL = 1, 2 scatterings and thermal effects [2]. Including the 
contributions of AL = 1 scatterings both in the wash-out term and in the CP asymmetry 
does not change our results significantly. We have summed up the contributions of the 
two quasi-degenerate RH neutrinos employing the property that eAr-^ = catj = e. This 
follows from having assumed that the decay rates of the two RH neutrinos are nearly 
equal and therefore both CP asymmetries are resonantly enhanced. If the decay rates 
are significantly different, then one should pick up only the CP asymmetry contribution 
which is resonantly enhanced. 

Finally, the time-dependent CP asymmetry relevant for resonant leptogenesis is given 
by US] 



2 / Kz' AM Tn, \ Tn. . Kz' AM V 
2sm — — -- — — — — —sm ' 



Im 



4 Fa/, / AM \ 2 F 



(AAt)^^ (AA% 1 + (AM/F^,; 
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(3) 



The CP asymmetry therefore consists of two blocks. The first one is the constant piece 
which is the usually adopted CP asymmetry and is resonantly enhanced for AM = Fat^. 
It allows efficient generation of the lepton asymmetry even for Mi ~ M2 as low as the 



TeV scale. The other block is made of two oscillating functions. The typical timescale for 
the variation of the CP asymmetry is 

1 ^2 ^^2 1 

t = = = ~ . (4) 

2H 2H{Mi) 2rjvi AM ^ ' 

The CP asymmetry grows for t < 1/AM and manifests its oscillation pattern only for 
t > 1/AM. The reader familiar with CP violation in neutral meson systems would 
prompltly recognize that the oscillation pattern originates from the CP violating decays 
of the two mixed states A^^i and A^2- These states do not propagate freely in the plasma 
though. If the timescale for the processes relevant for leptogenesis is much larger than 
~ 1/AM, the CP asymmetry should average to the constant value e quoted in the 
literature [20l [HI [21]. However, if the timescale of the evolution of the CP asymmetry is 
larger than or of the order of the timescale of the other processes, the time dependence of 
the CP asymmetry may not be neglected. This is precisely what happens in the resonant 
leptogenesis scenario where AM ~ Fatj ~ F^Vi. 

Since the strength of the interaction rates is dictated by the parameter K, we expect 
that in the strong wash-out regime, K ^ 1, the effect of the time dependence of the 
CP asymmetry is negligible. Due to the rapidly oscillating CP asymmetry, the lepton 
asymmetry should also rapidly oscillate and - at large times - reproduce the value usually 
quoted in the literature. On the contrary, in the case of weak or mild wash-out, K <1, 
the effect of the time dependence of the CP asymmetry should be magnified since the CP 
asymmetry oscillates with a period comparable to the time scale of the other interactions. 

The numerical solutions of the Boltzmann equations support these expectations. 
Figs. [T] and [2] show the evolution of the lepton asymmetry with and without the time 
dependence in the CP asymmetry for two representative cases of strong and weak wash- 
out, respectively. Fig. [3] shows the final baryon asymmetry computed taking into account 
and neglecting the time dependence in the CP asymmetry as a function of K. 

The numerical results can be analytically reproduced. In the strong wash-out regime, 
K ^ 1, {Yn — Y^) ~ {zK2{z) / AKg^) and the lepton asymmetry reads 

1 r'^ 



Y^^— dze(z)z'Ki(2)e-tJ. d.'(^'rKi(.')_ (5) 

Using the stepeest descent method, one can easily show that the final lepton asymmetry 
is equal to the one computed neglecting the time dependence in the CP asymmetry up to 
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Figure 1. The absolute value of the lepton asymmetry with the time dependence in 
the CP asymmetry (blue) and without it (green), as a function of z, for AM/Tn2 = 
TnJTn, = 1 and for K ^ 10. 




z 



Figure 2. The absolute value of the lepton asymmetry with the time dependence in 
the CP asymmetry (blue) and without it (green), as a function of z, for AAZ/FjVa — 
TnJTn^ = 1 and for K = lQ-\ 



small corrections 

The weak wash-out regime is much more interesting. Let us first remind the reader 
what happens in the usual case where the time dependence of the CP asymmetry is 
neglected. The final lepton asymmetry results from a cancellation between the (anti-) 
asymmetry generated when RH neutrinos are initially produced and the lepton asymmetry 
produced when they finally decay. It is useful to define the value Zec^ ^ 1 as the 
"time" when the A^-abundance reaches the equilibrium abundance: Yn^Zcc^) = Y^{zcq). 
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Figure 3. The absolute value of the final lepton asymmetry with the time dependence 
in the CP asymmetry (blue) and without (green), as a function of K, for AM/Tn^ = 



Since J^"'^ dz'{z')^Ki{z') ~ 37r/2, one finds that z^q is defined implicitly by the relation 
2^eq^e~^°i ^ {3'n-K/2). For z < Zeq, inverse decays dominate over decays and <C Y^^. 
From Eq. ([2]) one finds 

For z > Zeq, decays dominate over inverse decays, Y^ ^ Y^, and 

Y, F,(ze,)e--^i.'^^'^'^^(^')/^^(^') ^ Y,izeq)e~-/^^^''<\ (8) 
The lepton asymmetry is therefore given by 



2eK 



^9 



dz'(;^')'Ki(z')e-^^^""^"(^")'^^(^") 



37tK 
8g* 



-K r Az"z"Kx(z")/K2(z") 



e — e ~ - 1 + - — ~ . (9) 

Let us now come back to the case in which the time dependence of the CP asymmetry 
is accounted for. This time the near-cancellation between the (anti-) asymmetry generated 
when RH neutrinos are initially produced and the lepton asymmetry produced when they 
finally decay is not expected to hold since the asymmetry is modulated by the varying CP 
asymmetry. Again, we split the final lepton asymmetry as the sum of two contributions 
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2 4 6 8 10 

AM/rN2 

Figure 4. The absolute value of the final lepton asymmetry as a function of /S.M /Tm2 
for K — 10^^, Vm^I^Ni = 1- The blue solid line is obtained by numerically integrating 
the Boltzmann equations; the black dashed line represents the analytical approximation 
given in the text. The normalization of the asymmetry by (e/<?*) has been performed 
only with the part of e which is independent of /S.M /Tn2i so that the full dependence of 
Yc on AM/FjVj is explicitly shown. 
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(10) 



where in the first contribution we have approximated the CP asymmetry by e{z) ^ 
—e{Kz'^/2). This formula reproduces fairly well the complicated pattern shown in Fig. H] 
and shows that the resonance is displaced from the position AM/T^^ = 1 obtained 
when the time dependence of the CP asymmetry is neglected. Similarly to Eq. iQ, the 
lepton asymmetry scales like K^. However, this scaling is not due to cancellations among 
the asymmetries at different stages, but rather to the fact that for tiny values of K, or 
equivalently for t < 1/AM, the CP asymmetry is suppressed. 

For AM/r^Va = 1, Eq. ^ reduces to 



Yr ~ —K 



K H sm^ 

8 2 



Kz^ 




Figure 5. The ratio of the absolute value of the final lepton asymmetry with and 
without the time dependence in the CP asymmetry as a function of for AAf/FAr^ = 



For example, for K = 10~^, the ratio between the final lepton asymmetries given in 
Eqs. (fTTi) and ([9]) is ~ 11, which is in good agreement with our numerical results, as 
shown in Fig. [5l 



3. Resonant leptogenesis revisited: the flavoured case 

Let us now go beyond the one-fiavour approximation. As we have already mentioned, 
it is rigorously correct only when the interactions mediated by charged lepton Yukawa 
couplings are out of equilibrium. Supposing that leptogenesis takes place at temperatures 
T ~ Ml ~ M2, the one-fiavour approximation only holds for Mi ~ M2 > 10^^ GeV. In 
this range all the interactions mediated by the charged lepton Yukawa couplings are out 
of equilibrium and there is no notion of fiavour. One is allowed to perform a rotation in 
fiavour space to store all the lepton asymmetry in one fiavour, the total lepton number. 
However, at T ~ 10^^ GeV, the interactions mediated by the charged tau Yukawa come 
into equilibrium followed by those mediated by the charged muon Yukawa at T ~ 10^ 
GeV and the notion of fiavour becomes physical. Including the issue of fiavour can 
significantly affect the result for the final baryon asymmetry [22| [23t [2il[25l fT0l[26| [27] . 
Thermal leptogenesis is a dynamical process, involving the production and destruction 
of RH neutrinos and of the lepton asymmetry that is distributed among distinguishable 
fiavours. The processes which wash out lepton number are fiavour dependent, e.g. the 
inverse decays from electrons can destroy the lepton asymmetry carried by, and only by. 
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the electrons. When flavour is accounted for, the final value of the baryon asymmetry 
is the sum of three contributions. Each term is given by the CP asymmetry in a given 
flavour properly weighted by a wash-out factor induced by the lepton number violating 
processes for that flavour. 

The Boltzmann equations with flavour taken into account are 



1 



Yl = -2e,{z)Y^-^K,z'K,{z)Ye,, (^ = e,/x,r), (12) 



where 



2sm — — — — — — — — sm ' 



4 Tn, TnJ am V 2 r^, Tn, 

~ (AAt),,(AAt),, l + {AM/r^,r 

Notice that, for simplicity, we have again assumed that the partial rates for the decay of 
the two quasi-degenerate RH neutrinos are nearly equal. Extending our study to the case 
where the rates are different is straightforward. 

What we have learned so far is that significant differences with respect to the case 
in which the CP asymmetries are constant in time are expected in the weak wash-out 
regime. Let us then suppose that all flavours are in the weak wash-out regime and also 
that K < 1. If the time dependence of the CP asymmetry is neglected each individual 
flavour asymmetry suffers the usual cancellation and is given by Y^. ^ 2.8{ei/ g^)KKi [26] . 
Solving the Boltzmann equations along the same lines leading to the results (fTOj) and (fTTi) 
shows that flavour asymmetries, when the time dependence in ej is taken into account, 
are given by (for AM/T^^ = 1) 



Y. ~ 



457r Svr . 2 / ^< 

K H sm ' 

8 2 



(14) 



9* 

Therefore, the flavour asymmetries are increased by a factor ~ 10{K/Ki) compared to 

I We neglect both quantum flavour correlations [10 and the corrections arising from connecting the 
asymmetries in the lepton doublets to the asymmetries in the charges — B /S — Li conserved by weak 
sphalerons [22l HOl [251 126] . They introduce small corrections to our results. 
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Figure 6. The ratio of tlie absolute value of the final lepton asymmetry in the e flavour, 
with and without the time dependence in the CP asymmetry, as a function of Kf., for 
— 0.05, AAf/r^v^ = ^N-il^Ni — 1- The red solid line is obtained by numerically 
integrating the Boltzmann equations; the black dashed line represents the analytical 
approximation given in the text. 
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Figure 7. The ratio of the absolute value of the final lepton asymmetry in the e flavour, 
with and without the time dependence in the CP asymmetry, as a function of ifg, for 
= 10, AM/r^Va = ^N^/^Ni — 1- The red solid line is obtained by numerically 
integrating the Boltzmann equations; the black dashed line represents the analytical 
approximation given in the text. 



what is obtained when the CP asymmetry is assumed to be a constant. This result is 
confirmed by the numerical study shown in Fig. O 

Let us now consider the case in which some flavour is in the strong wash-out regime, 
so that K ^ 1, but some flavour i is weakly coupled, Ki < 1. Notice that this case is 
peculiar to the scenario where flavour is taken into account; it has no analogue in the 
one-flavour regime. Now the lepton asymmetry stored in the flavour i without accounting 
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for the time-dependence of the CP asymmetry reads 





= e,K, / dzY^'^z'K.iz) ^ 0.8^ (15) 
Jo 9* 

where we have made use of the fact that, for ^ 1, Fjv ~ Y^. The result is proportional 

to the weak wash-out parameter Ki because, in the limit of vanishing Ki, the asymmetry 

produced by inverse decays is cancelled by the one generated by the decays of the RH 

neutrinos. If the time dependence of the CP asymmetry is accounted for, the lepton 

asymmetry in the flavour i is given by 

Y.c^ -2 / dze,(z)F^e-^/^^'^^'(^')'^^(^') 







noo POO 

-2/ dze[iz)Y^ + Kj dze,iz)Y^z^K^{z). (16) 
Jo Jo 

where in the second line we have neglected the damping factor. The piece proportional 
to the derivative of the CP asymmetry may be evaluated as follows. Since K ^ 1, the 
RH abundance can be approximated as 

^fl-e-f^^^l (z<z,,. 



where Zeq — (6//^)^/^. The integral can now be evaluated as the sum of three different 
pieces, for < 2; < Zeq, Zeq < z < 1 and for 2; > 1. The latter is negligible since for large 
K the oscillating functions in the derivative of the CP asymmetry average to zero. The 
other two integrals can be easily evaluated using the expression ( ITTl) and give 



2 / dzeAzjYM ~ Zea-. -; cos ^ — sm ^ 

Jo 2(?, '""AM [AM \ 2 T^J \ 2 T^J 

The final lepton asymmetry in the flavour i is therefore 



'1^ 



Y. ^ O.S^K, - ^z,q^ 



2g, ^^AM 



r^, (Kzl^AM\ . (Kzl AM 
cos ^— — — sm ' 



AM \ 2 Tn.J V 2 



(19) 



For instance, for i^' = 10 and AM ~ F7V2, one obtains z^^ — 0.85 and the ratio between 
the asymmetries (|T9l) and (|T5l) goes like 1 + (0.19//^^), which agrees with the numerical 
results, as shown in Fig. [3 
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4. Conclusions 

Resonant leptogenesis has received much attention since it allows an efficient generation 
of the baryon asymmetry for RH neutrinos as light as the TeV scale. Through a non- 
equilibrium quantum field theory approach to leptogenesis, we have recently shown that 
the CP asymmetry parameter is not constant in time, but it varies with a typical timescale 
equal to the mass difference of the RH neutrinos [16]. In resonant leptogenesis, the two RH 
neutrinos Ni and N2 are almost degenerate in mass and the CP asymmetry from the decay 
of the ffist RH neutrino A^^i is resonantly enhanced if the mass difference AM = {M2 — M1) 
is of the order of the decay rate of the RH neutrinos. Therefore, the typical timescale 
of variation of the CP asymmetry can be larger than the timescale for the change of the 
abundances of the RH neutrinos. 

We have studied what differences arise with respect to the case in which the CP 
asymmetry is assumed to be a constant. Let us summarize our results: 

1) One-flavour case, valid for Mi ~ M2 > 10^^ GeV: the expression for the final 
baryon asymmetry differs from the results appeared so far in the litarature in the weak 
and mild wash-out regime, K < 1. The baryon asymmetry scales like K^, see Eqs. ( 1101) 
and (fTTj) . and is a factor (9(10) larger than the baryon asymmetry computed with constant 
CP asymmetry. 

2) Flavoured case, to be applied when Mi ~ M2 < 10^^ GeV: the expression for 
the final baryon asymmetry differs from the results appeared so far in the litarature for 
those flavours i which are in the weak and mild wash-out regime. If if < 1 and ifj < 1, 
the asymmetry in the flavour i scales like K"^, see Eq. f lT^ . and is enhanced by a factor 
~ 10{K/Ki) with respect to the case in which the CP asymmetry is constant. If K > 1 
and Ki < 1, the asymmetry in the flavour i is given by Eq. ([T^ and can be enhanced 
by a factor proportional to 1/Ki with respect to the case in which the CP asymmetry is 
constant. 

We conclude that the memory effects encoded in the time-dependent CP asymmetry 
may play an important role in resonant leptogenesis. 
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